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Introduction 

In  casting  of  single-crystal  turbine  blades  for  jet  engines,  the  formation  of  solidification  defects  has  become  an 
increasingly  important  problem  due  to  the  rising  levels  of  refractory  elements  in  Ni-based  superalloys.  Refractory 
elements,  which  are  beneficial  for  high-temperature  mechanical  properties,  enhance  density-driven  convective 
instabilities  underlying  the  formation  of  freckle  defects  in  directional  solidification.  In  support  of  an  effort  aimed 
at  the  development  of  validated  mathematical  criteria  for  predicting  solidification-defect  formation  in  superalloys, 
ab-initio  molecular  dynamics  simulations  have  been  performed  to  calculate  atomic  volumes  of  Ni-Al,  Ni-W,  Ni-Re 
and  Ni-Al-W  melts. 

Objective 

In  directional  solidification,  the  mass-density  difference  between  the  hot  liquid  near  the  solid-liquid  interface  and 
the  cooler  melt  at  the  top  and  side-walls  of  the  casting,  is  the  driving  force  for  convection  and  associated  formation 
of  solidification  defects.  This  mass-density  difference  originates  from  both  the  composition  and  temperature 
dependencies  of  the  liquid-phase  molar  volume  (R(c,  T)).  The  objective  of  this  work  is  the  application  of  ab-initio 
molecular  dynamics  simulations  as  a  computational  framework  for  aiding  the  development  and  validation  of  accurate 
models  for  V{c,T)  in  Ni-based  superalloys. 

Methodology 

The  properties  of  molten  Ni-based  alloys  have  been  computed  within  the  framework  of  ab-initio  molecular  dynamics 
(AIMD)  simulations.  For  500-atom  simulations,  ionic  positions  were  evolved  in  time  employing  classical  molecular 
dynamics  in  a  constant-temperature  and  volume  ensemble,  using  interatomic  forces  computed  directly  within  the 
framework  of  electronic  density  functional  theory  (DFT).  This  implementation  of  AIMD  made  use  of  the  commercial 
DFT  software  VASP  (Vienna  ab-initio  simulation  package).  Nearly  linear  scaling  is  achieved  for  500-atom  calculations 
run  on  up  to  approximately  128  processors  on  an  SGI  03K  and  Altix. 

Results 

Simulations  were  performed  for  elemental,  binary  and  ternary  alloys  of  Ni  with  Al,  W  and  Re  to  compute  equations  of 
state  at  temperatures  of  1830  and  1750  K.  Where  comparisons  with  measurements  are  available,  AIMD-calculated  vol¬ 
umes  agree  to  within  0.5-2. 5  %  of  experiment.  The  results  are  used  to  test  the  accuracy  of  the  predictions  of  a  recently 
proposed  parametrization  for  composition  and  temperature  dependent  molar  volumes  in  liquid  multicomponent  su¬ 
peralloys.  For  Ni-Re  the  model  is  found  to  be  highly  accurate  over  a  wide  range  of  compositions  while  for  Ni-W  it  gives 
rise  to  qualitatively  incorrect  predictions  at  higher  W  concentrations.  The  AIMD  results  thus  establish  the  limits  of 
the  accuracy  of  the  types  of  density  models  used  in  commercial  casting  software  and  are  suggest  ways  to  improve  them. 

Significance  to  DoD 

All  weapons  systems  which  employ  Ni-based  super-alloys  in  propulsion,  energy  generation,  temperature  protection 
systems  or  other  high  temperature  applications  will  benefit  from  this  research.  Relevant  systems  include:  current 
and  future  turbine  powered  aircraft  and  space  lift /re-entry  systems.  The  problem  of  freckle  formation  is  particularly 
important  for  the  development  of  superalloys  utilized  for  turbine-blades  in  high  performance  jet  engines. 
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In  the  casting  of  single-crystal  turbine  blades,  the  composition  (c)  and  temperature  (T)  dependen¬ 
cies  of  the  liquid-phase  molar  volume  (V{c,T))  play  a  critical  role  in  driving  convective  instabilities 
and  the  associated  formation  of  solidifcation  defects.  In  support  of  a  current  effort  aimed  at  the 
development  of  validated  mathematical  criteria  for  predicting  the  formation  of  solidification  defects 
in  Ni-based  superalloys,  ab-initio  molecular  dynamics  (AIMD)  simulations  have  been  performed 
for  elemental,  binary  and  ternary  alloys  of  Ni  with  Al,  W  and  Re,  to  compute  equations  of  state 
at  temperatures  of  1830  and  1750  K.  Where  comparisons  with  measurements  are  available,  AIMD- 
calculated  volumes  agree  to  within  0.5-2. 5  %  of  experiment.  The  results  are  used  to  test  the  accuracy 
of  the  predictions  of  a  recently  proposed  parametrization  for  composition  and  temperature  depen¬ 
dent  molar  volumes  in  liquid  multicomponent  superalloys.  For  Ni-Re  the  model  is  found  to  be  highly 
accurate  over  a  wide  range  of  compositions  while  for  Ni-W  it  gives  rise  to  qualitatively  incorrect 
predictions  at  higher  W  concentrations. 


I.  INTRODUCTION 

In  the  casting  of  single-crystal  turbine  blades  for  high-performance  jet  engines,  the  formation  of  solidification  defects 
has  become  an  increasingly  important  problem  due  to  the  rising  levels  of  refractory  elements  in  commercial  Ni-based 
superalloys.  While  refractory  elements  have  the  desirable  effect  of  improving  high-temperature  creep  strength,  they 
are  also  known  to  enhance  density- driven  convective  instabilities  underlying  the  formation  of  so-called  freckle  defects 
(chains  of  small  equiaxed  grains)  in  directional  solidification.  The  high-angle  grain  boundaries  and  composition 
gradients  associated  with  freckles  adversely  affect  mechanical  properties,  and  the  optimization  of  superalloys  for  high- 
performance  applications  thus  requires  quantitative  models  for  predicting  the  conditions  leading  to  their  formation. 

Freckles  develop  due  to  thermo-solutal  convection  in  the  molten  alloy  during  solidification,  and  it  has  been  suggested 
that  they  will  form  when  a  critical  Rayleigh  number  is  exceeded^  The  Rayleigh  number  (R)  is  a  measure  of  the 
ratio  of  the  buoyancy  force  to  the  retarding  frictional  force  in  the  mushy  zone: 

R  =  {Ap/p)gKl/aiy  (1) 

where  I  is  an  appropriate  length  scale,  K  is  the  average  permeability,  g  gravitational  acceleration,  a  is  the  thermal 
diffusivity,  u  is  the  kinematic  viscosity  of  the  fluid,  and  (Ap/p),  the  density  contrast,  is  a  measure  of  the  density 
variation  over  the  mushy  zone  between  the  solid  and  liquid  phases.  The  tendency  towards  defect  formation  is  highly 
sensitive  to  the  primary  dendrite  arm  spacing  (which  determines  the  permeability),  as  well  as  the  density  contrast. 
An  open  question  concerning  the  predictive  power  of  the  Rayleigh- number  criterion  for  freckle  formation  is  the  fact 
that  it  assumes  the  density  gradients  to  be  aligned  with  gravity.  This  is  not  the  case  in  most  castings  since  heat  losses 
through  the  mold  create  horizontal  thermal  gradients  and  hence  tilted  solid-liquidvinterfaces.  The  work  presented  in 
this  paper  is  a  component  of  a  multi-instituional  effort,  supported  through  the  Air-Force- Office  of  Scientific  Research 
Materials  Engineering  for  Affordable  New  Systems  (AFOSR-MEANS)  program,  aimed  at  the  development  of  validated 
mathematical  criteria  for  predicting  solidification-defect  formation  in  the  casting  of  single-crystal  Ni-based  superalloys. 
A  model  ternary  Ni-Al-W  superalloy  system  has  been  chosen  as  an  initial  focus  of  the  experimental  and  modeling 
efforts  to  validate  the  physical  models. 

In  the  validation  of  mathematical  models  for  freckle  formation  in  specific  alloy  systems,  it  is  critical  that  the  values  of 
the  various  parameters  in  Eq.  1  are  accurately  known.  Determination  of  the  permeability  and  density  gradients  poses 
a  major  challenge  in  this  regard.  In  the  MEANS  project,  values  for  the  permeability  are  being  derived  from  fluid- flow 
simulations  employing  realistic  models  for  the  mushy-zone  topology  obtained  from  three-dimensional  reconstructions 
derived  by  serial  sectioning^.  The  focus  of  the  work  presented  here,  performed  with  HPC  Challenge  support,  is 
the  application  of  ab-initio  molecular  dynamics  (AIMD)  simulations  as  a  computational  framework  for  aiding  the 
development  and  validation  of  accurate  models  for  Ap.  In  directional  solidification,  this  parameter  represents  the 
mass-density  difference  between  the  hot  liquid  near  the  dendritic  solid-liquid  interface  and  the  cooler  melt  at  the 
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top  and  side- walls  of  the  casting.  The  density  contrast  is  therefore  related  both  to  the  temperature  gradients  present 
during  casting,  as  well  as  the  variations  in  solute  concentration  in  the  melt  that  develop  as  a  consequence  of  equilibrium 
partitioning  between  solid  and  liquid  phases.  The  mass-density  difference  thus  originates  from  both  the  composition 
(c)  and  temperature  (T)  dependencies  of  the  liquid-phase  molar  volume  (^(c,  T)). 

Very  recently,  Mukai  et  have  performed  detailed  measurements  of  the  densities  of  several  binary  liquid  Ni- 

based  alloys,  and  a  few  representative  ternaries,  as  functions  of  c  and  T.  The  results  of  their  measurements  have  been 
used  to  develop  the  most  thorough  parameterization  for  V(c,  T)^  in  superalloys.  Due  to  the  lack  of  extensive  data 
for  multicomponent  systems,  a  key  assumption  made  in  the  developing  the  parameterization  is  that  the  partial- molar 
volume  of  each  solute  species  is  independent  of  the  compositions  of  the  others.  The  effects  upon  V(c,  T)  arising  from 
the  interactions  between  solute  atoms  of  different  types  are  thus  neglected.  Further,  measured  density  data  for  some 
key  elements  such  as  Re  are  unavailable  and  the  parameters  for  these  elements  are  treated  as  fitting  parameters  to  fit 
available  data  for  multicomponent  superalloys.  Comparisons  between  the  predictions  of  the  Mukai  parameterization 
and  available  data  for  densities  of  commercial  superalloys  shows  agreement  at  the  level  of  a  few  percent.  At  present 
experimental  data  for  ternary  Ni-Al-W  alloys  is  not  available,  and  estimates  of  the  accuracy  of  the  Mukai  model  for 
liquid  density  in  this  system,  which  forms  the  focus  of  our  MEANS  project,  are  thus  not  readily  obtained.  The  work 
presented  here  makes  use  of  state-of-the-art  AIMD  simulations  as  a  means  for  testing  the  accuracy  of  the  model  for 
parametrizing  p(c,  T)  in  this  system.  These  calculations  are  also  being  used  to  independently  test  the  accuracy  of  the 
Mukai  parameters  for  elements  where  direct  measurements  of  liquid  density  are  unavailable. 


II.  METHOD 

The  properties  of  molten  Ni-based  alloys  have  been  computed  within  the  framework  of  ab-initio  molecular  dynam¬ 
ics  (AIMD)  simulations.  In  these  simulations,  ionic  positions  were  evolved  in  time  employing  classical  Nose-Hoover 
dynamics  (fixed  particle  number,  volume  and  temperature  MD)^’^^  using  interatomic  forces  computed  directly  within 
the  framework  of  electronic  density  functional  theory  (DFT).  The  current  implementation  of  AIMD  makes  use  of 
the  commercial  DFT  software  VASP  (Vienna  ab-initio  simulation  package)  developed  at  the  Institut  fiir  Mate- 
rialphysik  of  the  Univeritat  Wien.  The  VASP  code  scales  well  for  systems  with  large  numbers  of  atoms  on  parallel 
computers  with  high-bandwidth  communications  networks.  The  computational  problem  involves  repeated  applica¬ 
tions  of  fast-fourier-transform  and  iterative  diagonalization  algorithms.  While  the  problem  is  far  from  embarrassingly 
parallel  overall,  nearly  linear  scaling  is  achieved  for  representative  calculations  up  to  approximately  128  processors  for 
AIMD  runs  performed  on  an  SGI  03K  and  Altix,  as  we  have  documented  previously The  underlying  algorithms 
employed  in  VASP  have  been  extensively  tested  and  optimized  for  parallel  processing  with  large  numbers  of  processors. 
The  production  runs  performed  in  this  work  have  typically  involved  the  use  of  64  processors. 

In  the  liquid-alloy  simulations  we  employed  500-atom  cubic  supercells  and  periodic  boundary  conditions  (with  fixed 
periodic  lengths).  Below  we  present  results  obtained  to  date  for  seven  alloy  compositions,  which  are  referred  to  in 
terms  of  the  number  of  each  type  of  atom  in  the  simulation  cell:  Nisoo,  Ni4ooAlioo,  Ni473W27,  Ni4ooWioo,  Ni473Re27, 
NUooReioo,  and  Ni436Al5oWi4.  Consistent  with  the  goals  summarized  in  the  previous  section,  the  alloy  compositions 
considered  in  this  work  are  chosen  to  test  the  accuracy  of  the  available  parameterization  for  V(c,  T)^  in  its  application 
to  the  ternary  Ni-Al-W  alloys  being  considered  for  the  validation  experiments  within  our  MEANS  program,  and  for 
important  superalloy  refractory  additions  such  as  Re,  for  which  experimentally  measured  data  is  currently  unavailable. 

The  details  of  the  calculations  are  as  described  by  us  previously^^  and  are  only  briefly  summarized  here.  Nose- 
Hoover  equations  of  motion  were  integrated  specifying  thermostat  temperatures  of  T=1830  and  T=1750  K.  The  time 
steps  (At)  for  the  MD  simulations  were  chosen  as  At=0.002  ps  or  At=0.003  ps,  to  ensure  that  the  conserved  energy  in 
the  Nose- Hoover  dynamics  displayed  a  drift  in  time  no  larger  than  1  meV/atom-ps.  The  starting  configuration  for  the 
Nisoo  and  NUooAlioo  simulations  were  obtained  from  snapshots  derived  from  classical  molecular-dynamics  simulations 
based  on  classical  interatomic  potentials'^;  for  the  Ni-W,  Ni-Re  and  Ni-Al-W  simulations  the  initial  configurations 
were  obtained  by  replacing  Ni  and  A1  atoms,  selected  at  random,  by  W  and  Re.  At  each  MD  time  step  the  self- 
consistent  charge  density  and  associated  interatomic  forces  were  computed  for  integration  of  the  ionic  trajectories,  and 
relevant  system  properties  including  the  pressure,  energy,  temperature  and  atomic  coordinates  were  stored.  The  DFT 
calculations  made  use  of  ultrasoft  pseudopotentials^^’^^  and  the  PW91  generalized-gradient  approximation  (GGA)^^. 
The  electronic  wavefunctions  were  represented  in  a  plane- wave  basis  set  with  a  cutoff  energy  of  260  eV.  The  simulations 
employ  a  single  k-point  (T)  in  performing  reciprocal-space  summations,  with  electronic  eigenvalues  occupied  according 
to  Fermi  statistics  with  an  electronic  temperature  equal  to  that  for  the  ions.  For  each  alloy  composition  simulations 
were  performed  over  a  range  of  volumes  spanning  plus  or  minus  five  percent  of  the  experimentally  measured  or 
estimated  equilibrium  value.  Total  simulation  times  ranged  from  5-10  ps,  a  time  interval  long  enough  to  obtain  good 
statistical  precision  in  calculated  equation-of-state  parameters'^.  As  described  below,  the  system  size  and  simulation 
time  were  sufficient  to  observe  appreciable  solute  diffusion  and  realistic  liquid  structures. 
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FIG.  1:  Calculated  mean-square  displacements  versus  simulation  time  derived  from  AIMD  simulations  for  elemental  Ni  and 
binary  Ni-X  (X=A1,  W  and  Re)  alloys  at  T=1830  K.  The  filled  areas  show  the  deviation  in  the  “random  walk”  from  the 
long-time  linear  diffusion  relation  (\r{t)  —  r(0)|^)  =  6Dt  for  each  species.  Estimates  of  the  self-diffusion  constants  (D),  reported 
in  units  of  jps  (10“^cm^/s),  are  indicated  in  each  of  the  panels. 


III.  RESULTS  AND  ANALYSIS 

The  main  focus  of  the  AIMD  simulations  is  the  calculation  of  molar  volumes  for  representative  binary  and  ternary 
alloy  compositions  in  Ni-based  superalloys.  In  calculations  of  such  equilibrium  properties  for  a  liquid  alloy,  it  is  impor¬ 
tant  to  establish  that  (i)  the  simulations  are  run  sufficiently  long  in  time  to  observe  appreciable  solute  interdiffusion, 
and  (ii)  that  the  periodic  system  size  is  large  enough  to  give  a  realistic  liquid  structure.  If  the  runs  are  performed 
for  too  short  a  period,  the  calculated  thermal  averages  may  be  biased  by  the  choice  of  the  initial  atomic  configura¬ 
tions.  Further,  if  the  system  is  too  small  the  positional  correlations  may  not  decay  sufficiently  over  the  length  of  the 
simulation  box,  giving  rise  to  an  artifically  enhanced  degree  of  order  in  the  simulated  melt  structures. 

Considering  first  point  (i).  Fig.  1  plots  averaged  mean-square  displacements  for  Ni,  W,  A1  and  Re,  as  calculated 
from  the  T=1830  K  simulations  for  the  pure-Ni  melt,  and  for  each  of  the  binary  NUoqXioo  samples  (X=A1,  W,  and 
Re).  The  slopes  of  the  best-fit  lines  relating  mean-square  displacements  and  time  are  used  to  derive  the  estimates 
of  the  concentration-dependent  tracer-diffusion  constants  listed  in  the  panels.  A  comparison  of  the  results  indicates 
that  Re  and  W  are  the  slowest  diffusing  elements  in  these  alloys,  as  expected  based  on  their  large  sizes  and  heavy 
atomic  masses.  It  is  also  interesting  to  note  that  the  Ni  diffusion  constant  is  lower  in  the  Ni-X  melts,  relative  to 
pure  Ni,  despite  the  fact  that  the  average  atomic  volume  is  larger  in  the  alloys;  this  result  may  be  a  reflection  of 
relatively  strong  solute- solvent  interactions  in  the  alloy  melts  (see  below).  From  the  standpoint  of  the  accuracy  of  the 
calculations,  it  is  important  to  note  that  Fig.  1  establishes  that  significant  interdiffusion  is  occuring  over  the  timescale 
of  the  AIMD  simulations.  Thus,  it  is  expected  that  the  runs  are  performed  over  sufficiently  long  times  to  give  robust 
estimates  of  the  equilibrium  molar  volumes.  This  expectation  is  supported  by  the  relatively  small  estimated  statistical 
uncertainties  in  this  property  quoted  below.  It  is  further  supported  by  the  results  of  a  test  performed  for  binary  Ni-W 
samples,  whereby  two  different  initial  atomic  configurations  were  used  to  compute  the  equilibrium  volume  and  the 
final  results  were  found  to  agree  within  statistical  uncertainties. 

Figure  2  plots  calculated  radial  distribution  functions  gQ,_/3(r)  for  each  of  the  pairs  of  species  in  pure  Ni 

and  binary  NUoqXioo  (X=A1,  W  and  Re)  melts.  Considering  point  (ii)  above,  it  is  noteworthy  that  each  of  these 
distribution  functions  converges  to  the  uncorrelated  limit  of  unity  over  a  distance  less  than  half  the  periodic  dimensions 
of  our  simulation  cells.  The  results  in  Fig.  2  also  give  insight  into  the  nature  of  the  chemical  ordering  in  the  alloy 
melts.  Specifically,  the  larger  value  of  the  main  peak  g^i-x,  relative  to  those  for  the  solute-solute  (X-X)  and  solvent- 
solvent  (Ni-Ni)  pairs,  is  indicative  of  a  preference  for  neighbors  between  unlike  species,  i.e.,  for  chemical  short-range 
ordering.  Similarly,  it  is  noteworthy  that  the  average  nearest- neighbor  bond  length  for  Ni-X  pairs  (approximately 
2.485  A)  is  significantly  smaller  than  the  average  of  the  Ni-Ni  and  X-X  bond  lengths  (roughly  2.59  A).  Both  results 
are  indicative  of  a  significant  preference  for  chemical  short-range  ordering  between  Ni  and  X  species  in  each  of  the 
systems  considered. 

To  compute  equilibrium  values  of  U,  AIMD  simulations  were  performed  for  a  series  of  volumes  at  each  composition 
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FIG.  2:  Calculated  radial  distribution  functions  derived  from  AIMD  simulations  for  elemental  Ni  and  binary  Ni-X  (X=A1,  W 
and  Re)  alloys  at  T=1830  K.  Results  are  plotted  solvent-solvent  (Ni-Ni),  solute-solven  (Ni-X)  and  solvent-solvent  (X-X)  pairs. 


TABLE  I:  Calculated  atomic  volumes  for  molten  Ni-based  alloys  at  T=1830  K  and  T=1750  K.  The  numbers  in  parentheses 
represent  estimated  95%  confidence  intervals  on  the  last  digit.  The  calculated  (AIMD)  results  are  compared  with  the  predictions 
of  the  parametrized  model  due  to  Mukai  et  al.^ . 


Composition 

V(T=1750  K)  (cmVmole) 

AIMD  Mukai 

V(T=1830  K) 

AIMD 

(cm^/mole) 

Mukai 

Nisoo 

7.64(1) 

7.4597 

7.71(1) 

7.5727 

NuooAlioo 

8.01(1) 

7.7628 

7.94(1) 

7.8525 

NuooWioo 

6.5192 

8.05(1) 

7.3598 

Ni473W27 

7.5534 

7.70(1) 

7.6557 

NRooRcioo 

7.5869 

7.98(1) 

7.8695 

Ni473Re27 

7.4940 

7.76(1) 

7.6557 

Ni436Al5oWl4 

7.6920 

7.81(1) 

7.7682 

and  the  average  pressure  (P)  was  calculated  as  the  basis  for  fitting  the  equation  of  state.  As  described  previously 
for  each  simulation  at  a  given  volume  the  results  show  that  after  an  intial  short  transient,  lasting  a  few  tenths  of  a  ps, 
the  pressure  displays  equilibrium  fluctuations  about  a  well  defined  average  value.  An  analysis  of  the  time  scale  and 
amplitude  of  these  equilibrium  fluctuations  is  used  to  derive  estimates  of  the  statistical  uncertainties  in  the  calculated 
time-averaged  pressure.  For  each  temperature  and  composition  a  total  of  three  to  four  volumes  were  considered.  The 
calculated  time-averaged  pressures  were  then  used  to  derive  an  equation  of  state,  assuming  a  quadratic  relationship 
between  P  and  V.  From  the  fitted  equation  of  state  we  computed  the  equilibrium  (zero-pressure)  volume  (R),  as  well 
as  the  bulk  modulus  and  its  pressure  derivative^^.  Finally,  statistical  uncertainties  were  estimated  based  on  the  error 
estimates  for  each  of  the  individual  average  pressures.  The  results  are  recorded  in  Table  I,  listing  calculated  equilibrium 
volumes  for  each  of  the  seven  liquid  Ni-based  alloy  compositions  considered  to  date.  Data  for  all  compositions  at 
T=1830  K  are  included,  and  initial  results  for  two  compositions  are  also  included  at  T=1750  K. 

In  Table  I,  values  from  the  most  recent  parametrization  of  molten  superalloy  densities,  due  to  Mukai  et  are 
listed  along  side  the  AIMD  results  for  each  of  the  alloy  compositions  considered.  For  elemental  Ni,  and  the  binary 
Ni-Al  and  Ni-W  systems,  the  Mukai  model  reproduces  well  the  directly  measured  densities  obtained  at  these  same 
compositions  and  temperatures^’^.  Thus,  for  the  Nisoo,  NUooAlioo  and  Ni473W27  systems,  a  comparison  between  the 
Mukai  and  AIMD  results  can  be  considered  as  representing  a  direct  comparison  between  experimental  measurements 
and  calculations.  The  largest  discrepancy  is  found  for  elemental  Ni,  where  the  agreement  between  measured  and 
calculated  values  of  V  is  still  very  good,  at  the  level  of  1.8  %  at  T=1830  K,  and  2.4  %  at  T=1750  K.  For  Ni-Al 
and  Ni-W,  the  agreement  between  experiment  and  theory  is  slightly  better,  1.1%  and  0.6  %,  respectively.  For  each 
of  these  systems  the  calculated  molar  volumes  are  larger  than  the  measurements,  a  trend  that  is  often  found  for 
solid-state  ab-initio  calculations  performed  within  the  generalized-gradient  approximation  to  DFT.  As  a  result  of  the 
comparisons  between  experiment  and  theory  for  Ni,  Ni-Al  and  Ni-W  systems,  it  is  reasonable  to  conclude  that  the 
GCA  introduces  a  positive  systematic  error  in  the  prediction  of  the  Ni-based  liquid  atomic  volumes,  on  the  order  of 
0.5-2. 5  %  for  the  alloy  compositions  and  temperatures  considered. 

As  described  in  the  Introduction,  the  primary  motivation  for  the  present  work  is  to  use  AIMD  simulations  as  a 
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FIG.  3:  Molar  volumes  versus  concentration  in  binary  Ni-W  and  Ni-Re  melts  at  T=1833  K.  The  solid  blue  and  red  circles 
represent  results  from  AIMD  simulations  and  the  parametrization  due  to  Mukai  et  a/.®,  respectively.  The  solid  lines  represent 
best  fits  to  the  data,  assuming  a  quadratic  and  linear  dependencies  of  volume  on  concentration  for  Ni-W  and  Ni-Re,  respectively. 
The  insets  in  each  panel  give  calculated  ratios  between  the  partial  molar  volumes  of  the  solutes  (at  infinite  dilution)  and  the 
molar  volume  of  the  pure  Ni  melt. 


framework  for  developing  and  validating  predictive  models  for  p(c,  T)  in  Ni-based  superalloys.  For  this  purpose,  we 
consider  further  the  most  current  parametrization  for  V{c^T)  proposed  by  Mukai  et  al.^.  We  consider  cases  where  there 
is  no  experimental  data  available,  and  use  the  AIMD  results  as  a  basis  for  testing  the  accuracy  of  the  Mukai  model 
in  predictions  for  compositions  where  it  has  not  been  fit.  We  consider  first  a  ternay  Ni-Al-W  melt  representative 
of  the  compositions  that  are  being  considered  in  the  experimental  effort  within  our  MEANS  program.  The  alloy 
Ni436Al5oWi4  listed  in  Table  I  has  a  composition  which  represents  an  average  of  those  for  the  binary  Ni4ooAlioo  and 
Ni473W27  melts  also  considered  in  this  work.  Both  the  AIMD  and  Mukai  results  predict  that  the  volume  of  the 
ternary  alloy  is  well  approximated  (to  within  0.2  %)  by  a  linear  interpolation  between  the  binary  alloy  values.  This 
result  serves  suggests  that  one  of  the  main  assumptions  in  the  Mukai  parametrization,  namely  that  the  partial  molar 
volumes  of  the  binary  alloys  can  be  used  to  interpolate  molar  volumes  for  multicomponent  systems,  is  reasonably 
accurate  in  the  Ni-Al-W  system. 

A  second  test  of  the  Mukai  model  is  made  based  on  the  results  plotted  in  the  right  panel  of  Fig.  3,  corresponding 
to  composition  dependent  molar  volumes  for  Ni-Re  alloys.  In  this  binary  system  no  measured  data  is  available  for 
liquid  volumes,  and  the  partial  molar  volume  of  Re  was  obtained  as  a  fitting  parameter  by  Mukai  et  al.  to  reproduce 
data  from  a  multicomponent  superalloy  system.  The  comparison  between  the  AIMD  and  Mukai  results  for  Ni-Re  are 
seen  to  be  in  very  good  agreement:  both  sets  of  results  show  a  highly  linear  dependence  of  V  on  Re  concentration. 
From  the  slope  of  the  fitted  lines  in  Fig.  3  the  partial  molar  volume  of  Re  (VRe)  in  molten  Ni  can  be  extracted.  The 
AIMD  simulations  and  Mukai  model  are  in  very  good  agreement,  predicting  values  for  Vrc  that  are  approximately 
20  %  larger  than  the  molar  volume  of  molten  Ni  at  T=1830  K  (see  insets  in  Fig.  3).  While  the  comparison  between 
AIMD  and  the  Mukai  model  is  thus  very  favorable  for  Ni-Re,  stong  discrepancies  are  found  for  the  Ni-W  system, 
as  indicated  by  the  results  in  the  left  panel  of  Fig.  3.  For  Ni-W  both  the  Mukai  model  and  the  AIMD  calculations 
predict  strongly  non-linear  dependencies  of  V  on  composition,  but  qualitatively  different  trends  are  found  for  the 
nature  of  this  non-linearity.  Specifically,  the  Mukai  model  predicts  that  the  addition  of  high  concentrations  of  W 
gives  rise  to  a  deerease  in  the  volume  of  molten  Ni-W  alloys  relative  to  pure  Ni,  in  sharp  contrast  with  the  AIMD 
calculations.  The  slopes  of  the  fitted  volume  versus  concentration  relations  in  the  pure-Ni  limit  can  be  used  to  derive 
the  partial  molar  volume  for  W  in  Ni  at  infinite  dilution  (E^),  and  for  this  quantity  relatively  large  discrepancies 
are  seen  again  between  the  AIMD  values  and  the  predictions  of  the  Mukai  model.  The  AIMD  results  thus  call  into 
question  the  accuracy  of  the  Mukai  parametrization  for  concentrated  Ni-W  melts,  and  suggest  that  the  reliability  of 
this  parameterization  for  superalloys  containing  W  is  limited  to  the  range  of  dilute  W  compositions  where  it  was  fit. 


IV.  SUMMARY 


AIMD  simulations  have  been  performed  to  compute  liquid-phase  atomic  volumes  for  seven  elemental,  binary  and 
ternary  alloy  compositions  of  molten  Ni  with  Al,  W  and  Re  additions  to  compute  atomic  volumes  at  temperatures 
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of  1830  and  1750  K.  Comparisons  between  AIMD  results  and  direct  measurements  for  Ni,  Ni-Al  and  Ni-W  yield 
agreement  at  the  level  of  0.5-2. 5  %,  with  the  calculations  being  systematically  larger  than  experiment.  Simulations 
were  performed  for  a  ternary  Ni-Al- W  composition  where  the  accuracy  of  the  most  recent  parametrization  for  molten 
superalloy  densities  was  affirmed.  The  results  also  support  the  the  accuracy  of  the  parametrization  for  Ni-Re,  but 
suggest  that  its  reliability  for  superalloys  containing  W  is  limited  to  the  range  of  dilute  W  compositions  where  it  was  fit. 
Future  work  will  involve  further  AIMD  calculations  to  more  thorougly  test,  and  refine  as  necessary,  the  accuracy  of  the 
current  published  parametrization  for  V{c,T)  over  range  of  compositions  and  T  relevant  for  superalloy  solidification. 
This  will  include  calculations  for  Ni-Ta,  as  well  as  ternary  Ni-Al-Re  and  Ni-Al- Ta  alloys  to  test  (i)  the  reliability  of  the 
main  assumptions  underlying  the  published  parametrization,  related  to  the  interpolation  of  multicomponent  volumes 
from  binary  data,  and  (ii)  the  origin  of  the  anamolously  large  temperature  dependencies  found  in  the  parametrized 
values  of  the  partial  molar  volumes  for  the  refractory  elements  W,  Ta  and  Re^.  Future  calculations  will  benefit 
from  the  completion,  over  the  past  fiscal  year,  of  an  extensive  effort  leading  to  the  implementation  of  fixed-pressure 
dynamics  into  the  VASP  code^^,  which  should  allow  calculation  of  equilibrium  atomic  volumes  in  about  one-third 
the  time  required  by  fitting  equations  of  state  to  fixed-volume  simulation  data.  Further,  we  will  explore  statistical 
reweighting  techniques  to  evaluate  pressure  (or  volume)  responses  arising  from  dilute  solute  additions  to  molten 
Ni-based  superalloys. 
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